Spin-polaron band structure and hole pockets in underdoped cuprates 
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We present a variational approach based on the string picture to analyze the internal structure 
and dispersion of spin polarons with different symmetries in an antiferromagnet. We then use this to 
discuss the properties of underdoped cuprate superconductor within the 'doped insulator' picture. 
The theory explains the remnant Fermi surface for the undoped compunds, as well as hole pockets, 
Fermi arcs, high energy pseudogap and the the mid-infrared band in doped materials. Destructive 
interference between the phases of a photohole near F and the internal phases of the Zhang Rice 
singlet combined with our theory moreover explains the 'waterfall' phenomenon. 
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I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tivity (HTSC) in doped antiferromagnetic (AF) insula- 
tors, the research on that phenomenon concentrates to a 
great extent on the properties of single-particle-like exci- 
tations in such systems. It is also obvious that the micro- 
scopic mechanism of HTSC should hinge on the interac- 
tion between quasiparticles. At half filling the Hubbard 
model, which is a generic model used to describe strongly 
correlated systems, is an insulator for large enough C//i[l| 
whereas for low electron density the model is expected to 
be a Fermi liquid with a Fermi surface volume in accor- 
dance with the Luttinger theorem for any ?7/i[l|. With 
increasing doping one might therefore expect a phase 
transition from a 'correlation dominated' phase near half- 
filling to a Fermi-liquid phase for low density. The key 
property of the correlation dominated phase thereby is 
the splitting of the physical electron into the two Hub- 
bard bands which correspond to Fermionic holes and dou- 
ble occupancies moving in a 'background' of singly occu- 
pied sites, whereby the electrons forming the background 
retain only their spin degrees of freedom. This - and 
not a half-filled Fermi surface - is the picture underlying 
all successful theories for the angle resolved photoemis- 
sion spectroscopy (ARPES) data 3| obtained in insulat- 
ing compounds ^,y,i6|,|7|,|i,i,|ia[iil,[i3]. The question 
then is to what extent these Hubbard bands can be doped 
before the two-band structure collapses, the spin back- 
ground 'melts' into a Fermi sea and the Luttinger Fermi 
surface is regained. The key property of this 'doped in- 
sulator' phase should be a Fermi surface with a volume 
proportional to the number of doped holes because these 
are the only mobile Fermions. 

ARPES on cuprate superconductors has produced a 
wealth of information [is'l although it has proven difficult 
to extract a consistent picture. From spectra taken near 
optimal doping it was concluded that ARPES shows evi- 
dence for a 'large' Fermi surface consistent with the Lut- 



tinger theorem and band calculations |l3| which would 
imply that the doped Hubbard bands never exist. Sur- 
prisingly enough, however, insulating compounds such 
as Sr2Cu02Cl2 seemed to show a very similar Fermi sur- 
face as well - which has been termed the 'remnant Fermi 
surface' 14j . ARPES spectra from insulators show a band 
which disperses towards lower binding energy and then 
rapidly looses weight just as if it would cross a Fermi 
level. The only possible explanation for this phenomenon 
is a strong and systematic variation of the spectral weight 
of the conduction band which drops to near zero abruptly 
at a line in k-space which roughly coincides with the non- 
interacting Fermi surface. One can give simple arguments 
why such a behaviour is to be expected |.15<]. Since there is 
no reason why such a strong variation of spectral weight 
and bands with almost no spectral weight should occur 
only in insulators one should be cautioned that Fermi 
surface maps for doped compounds may not show the 
full picture either. 

In the underdoped compounds the upper Hubbard 
band, while rapidly loosing spectral weight, still can be 
clearly resolved ARPES shows rather structureless 
spectra which are usually interpreted in terms of a 'high 
energy feature' and a 'leading edge shift' or, alternatively, 
a high energy and low energy pseudogap [l3|- In any case 
there is definitely no 'large' Fermi surface, instead Fermi 
surface maps show 'Fermi arcs'. Bearing in mind the rem- 
nant Fermi surface in the insulators suggests to interpret 
these arcs as being the inner part of a hole pocket cen- 
tered on (§, f ) with the part of the pocket facing (tt, tt) 
having too small spectral weight to be seen in ARPES. 
Assuming that the electronic structure in the underdoped 
compound can be described in simplest approximation as 
holes doped into the quasiparticle band of the insulator 
moreover would give an explanation for the high energy 
pseudogap - the dispersion of the 'high energy feature' 
- in that it simply reflects the hole dispersion in the 
insulator p^. Such a 'strict' rigid band picture would 
not explain the low energy pseudogap or the tempera- 
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ture and doping dependence of either high or low energy 
pseudogap. On the other hand the closing of the pseudo- 
gap with increasing doping can be seen already in cluster 
simulations ^r?! and it has been pointed out there that 
the mechanism may be an effective downward renormal- 
ization of the t' and t" terms due to decrease in the spin 
correlation length. Moreover, the low energy pseudogap 
by its definition in terms of the leading edge shift has 
no immediate connection with a dispersion and may be 
determined e.g. by T-dependent linewidths as discussed 
by Storey et al.^. 

The compound Ca2-a;Naa;Cu02Cl2 shows rather clear 
evidence for the doped insulator picture [l^ in that 
the dispersion in the doped case is virtually identi- 
cal to that in the undoped compound and even the 
part of the quasiparticle band facing (tt, tt) may have 
been observed. Very recently ARPES experiments seem 
to have provided direct evidence for hole pockets in 
Lai.48Ndo.4Sro.i2Cu04^], with the part of the pockets 
facing (tt, tt) indeed having small spectral weight - as one 
would have expected on the basis of the remnant Fermi 
surface. 

It had been noted early on [2l| that hole pockets with 
a volume proportional to the hole concentration would 
explain the scaling of the low temperature Hall constant 
with hole concentration in the underdoped materials (2^ . 
[23I I and that the apparent discrepancy between the 'small 
Fermi surface' suggested by transport measurements and 
the 'large Fermi surface' seen in ARPES may be due to 
a systematic variation of the quasiparticle weight al ong 
the hole pocket. [2^. While the recent discovery [25l. [26L 
of Shubnikov-deHaas oscillations in some underdoped 
compounds initially seemed to provide strong evidence 
for hole pockets the finding that the oscillations are likel; 
caused by electron-like rather than hole-like pockets 
has complicated matters again. 

Lastly, exact diagonalization studies of the t-J model 
provide clear evidence that the Fermi surface for hole 
dopings around 10% takes the form of hole pockets [29I 
|30| . Careful analysis of exact diagonalization results 
shows that the single particle the spectral function for the 
doped t-J model is quite consistent with rigid band fill- 
ing of the quasiparticle band seen at half-filling, provided 
one takes into account the formation of hole pairs 24 1 . 
By calculating the spectral function for dressed hole 
operators jSlj rather than the bare electron operators is 
can moreover be shown that the quasiparticles in the 
doped system have very nearly the same internal struc- 
ture as in the undoped onejs^l- 

Here we take the point of view that the underdoped 
regime in high-temperature superconductors precisely 
corresponds to the doped insulator phase. We show that 
many properties of the underdoped phase - the remnant 
Fermi surface, the Fermi arcs, the high-energy feature 
seen in ARPES, the mid infrared band seen in the opti- 
cal conductivity - find a simple and natural explanation 
in the dispersion and internal structure of the quasipar- 
ticles which correspond to holes heavily dressed by spin 



excitations. Since there are indications that the heavily 
overdoped phase is essentially a Fermi liquid this would 
imply that the phase transition from the correlation dom- 
inated phase to the Fermi liquid phase occurs at optimal 
doping. This would then be a quantum phase transition 
where none of the two phases has any kind of order - 
rather they differ in the topology and volume of their 
Fermi surfaces. An indication of this transition can in 
fact be seen in the dynamical spin and density correla- 
tion function obtained by exact diagonalization of small 
clusters. In the underdoped regime spin and density cor- 
relation functions are very different and the density cor- 
relation function takes the form of extended incoherent 
continually This form of the density correlation function 
can be explained quantitatively within the string picture 
for a single hole 34] . For doping levels higher than opti- 
mal, spin and density correlation function become simi- 
lar and can be explained well as particle-hole transitions 
across an essentially free-electron-like Fermi surface[35]. 

While a theory for such a transition would be highly 
desirable but very challenging the present paper has a 
more modest goal: we want to show that many features 
of undoped and underdoped cuprates can be explained 
by a very simple theory which assumes continuity with 
the insulator. The calculation will be performed in the 
framework of the t-J model [1^, [s^] extended by terms 
enabling hopping to second and third nearest neighbors 
with hopping integrals t' and t" respectively jsS, 39]. We 
use standard values t = 0.35eV, t' — —0.12 eV, t" — 
0.08 eV, and J = 0.14eV chosen so as to reproduce the 
measured Fermi surface of hole doped cuprates for high 
doping levels. 



II. CONSTRUCTION OF LOCALIZED BASIS 
STATES 



Since we want to study the doped insulator we con- 
sider the motion of a single hole in an antiferromagnet- 
ically ordered 'spin background'. All processes analyzed 
in the following actually require only short range antifer- 
romagnetic correlations - one may therefore expect that 
hole motion in a state with short range antiferromagnetic 
order but no long range order will involve very similar 
processes so that e.g. the internal structure of the quasi- 
particles and the dispersion relation of a hole should not 
change drastically. The construction of spin polaron (SP) 
states including the excited states w as p erformed in sev- 
eral earlier publications [iO, Sll, [H, lU [11] . In order to 
make this paper self-contained we will now briefly repeat 
that construction. For definiteness we will assume that a 
|-spin has been removed from the system and study the 
motion of the resulting hole. We denote the J,-sublattice 
by A 

To begin with we define Hq = Ht + Hjsing to be the 
sum of the nearest neighbor hopping cx t and the longi- 
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FIG. 1: (Color online) The mechanism of the string effect. 
Slanted crosses represent links with the contribution to the 
Ising part of the exchange energy higher by J/2 as compared 
to the Neel state, (b), a string state obtained by a single 
move of the hole created at site i. (c), (d), states obtained 
respectively by two and three consecutive moves. 



tudinal part of the Heisenberg exchange: 
Ht 



Hi 



Sing 
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Hf) is frequently referred to as the t — model. In a 
first step, we seek approximate eigenstates of Hq which 
are localized due to the string effect. 

The mechanism of the string effect is shown in FigH] 
By creating a hole in the Neel state - i.e. the ground 
state of Hjsing - at site i and acting repeatedly with 
the hopping term we generate a basis of string states 
\Vi} where Vi = {i,j,...,n) is shorthand for the sites 
i,j,...,n visited by the hole. All spins on these sites 
have been displaced by one lattice spacing and thus are 
inverted relative to the Neel order. 

All [Pi) are eigenstates of Hi sing- Taking the energy 
of the string with length i.e. the bare hole as the zero 
of energy and denoting by v the length of the string - i.e. 
the number of shifted spins - the eigenvalue is 



1) 



(1) 



This is exact for v < 2 and is true for 'most' longer 
strings as well. For simplicity we assume ([T|) to be true 
for any string. This implies that the hole is trapped in a 
linearly ascending potential and all eigenstates of Hq are 
localized. The main deviations from ([1]) occur for 'loops' 
as discussed by TrugmanQ which in fact lead to hole- 
propagation even in the t — Jz model. Such loops pose no 



FIG. 2: Dependence of the sign of a given path on the direc- 
tion of the first move for the d^2_j^2-wave SP, (b), pa;-wave 
SP, (c), and Py-wave SP, (d). For completeness the schematic 
representation for the s-wave SP has also been shown, (a). 



fundamental problem for the present formalism and can 
be dealt with by introducing the co ncep t of 'irreducible 
paths' as discussed in detail in Ref. [4(j |. 

Next we note that by acting with a point group oper- 
ation, which leaves the initial site i invariant, the string 
states \Vi) are transformed into one another. We can 
therefore define linear combinations of the [Pi) which 
transform like the basis states of the irreducible represen- 
tations of Ciy under these point group operations. Then 
we make the following ansatz for a localized eigenstate of 
Hq 



1^ 



(o,m)\ _ 



a. 



(2) 



where o G {s,Px,Py,dx2-y2^dxy ■ ■ ■} denotes the sym- 
metry or 'orbital character' of the state and m labels the 
excitation number for a given symmetry. In keeping with 

^ we moreover assume that each coefficient a^.'™'' can 

be factorized into a sign (f>^'^ - which plays the role of 
an 'angular wave function' - and a 'radial wave function' 



which depends only on the length of the string: 



(o,m) 



(o,m) 



(3) 



For an Ai (s-wave) state the sign (j)-p. is obviously uni- 
form for all paths. For an E (p-wave) or Bi ((ia;2_j^2-wave) 
state (f>-p- is determined by the direction of the first hop 
away from the site i as shown in Figure [21 The string 
of length i.e. the bare hole at site i is invariant under 
all point group operations and hence has nonvanishing 
weight only in the s-like state - this implies that E and 
Bi states are higher in energy than the Ai state because 
the are composed of strings with length > 1 and hence 
a minimum of three frustrated bonds. For the two re- 
maining representations A2 (g-wave) and B2 (dxy-wave) 
it can be shown that only strings with a minimum length 
of 2 have nonvanishing weight in the corresponding SP 
states - these states therefore are even higher in energy 
and we omit them. For each symmetry sector o we can 
now set up the Schrodinger equation 



(4) 
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thereby assuming ((T]) and solve for the eigenenergies 
E^o,m) g^^^ ^j^g coefficients a - this is explained in Ap- 
pendix [A] As expected for a linearly ascending potential, 
the a's are rapidly decreasing with increasing length of 
the string even for the physical parameter range t/ J « 3. 
It turns out the Schrodinger equation for the coefficients 
a for the dj.'2_y2 SP is identical to that for the p-like ones. 

To shorten the notation we call the coefficients a^f''^^ for 
the lowest s-like SP a^, and instead of those for the p-like 
SP, Q;[f'°^ and (ia;2-y2-like SP, aif''^'' we use a'^/\/2 and 
respectively. 
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III. EFFECTIVE MULTI-BAND MODEL FOR 
SPIN POLARONS 

So far we have found SP states ([2]) which form a set 
of approximate localized eigenstates of Hq at each of the 
sites of the sublattice A. Next we note that the remain- 
ing part Hi of the i JM - which comprises the transverse 
part of the Heisenberg exchange and the hopping terms 
DC t' , t" - has nonvanishing matrix elements between SP 
states centered on neighboring sites i and j. One im- 
portant mechanism leading to such a matrix element is 
the truncation of the string shown in Figure [31 By flip- 
ping the flrst two spins of the defect string the starting 
point of the string is shifted to a second or third nearest 
neighbor while the length of the string is reduced by two. 
Since the coefficients a^.'™' and ai^, of the initial and 

final string are known by solution of ^ and the strength 
of the spin-flip term is J/2 the corresponding matrix el- 
ement is easily evaluated. Similarly, the hopping terms 
cx t' ,t" allow for the hopping of the bare hole between 
the sites of one sublattice (see Figure In addition 
there is the 'loop hopping' [3| and actually a wide variety 
of additional processes, which are discussed in Appendix 

Assuming that the matrix elements are known we de- 
fine Fourier transforms 



|vI.W(k))^7|^|*f))e-''-^^- 



(5) 



where A = (o,m). Next we make the LCAO-like ansatz 
for a propagating single-hole state 



|$,(k))=^i;i')(k)|vI/W(k)) 



(6) 



which leads to a generalized eigenvalue problem of the 
form 

i/,//(k)vKk) = i?Kk)Oe//(k)vKk) (7) 

where the Hamilton and overlap matrices are given by 

i?(AA') - <5;,vi?(^)+r(^^')(k), (8) 
O(AA') = <5AA'+0(^^'Hk) (9) 
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FIG. 3: (Color online) (a)^(b): a string of length 2 starting 
at i is truncated to a string of length at j by the transverse 
part of the Heisenberg exchange. 

(c) — > (d) a string of length 3 is reduced to one of length 1. 
(e) (f) a bare hole at site i ('string of length 0') is trans- 
ported to j by t'-hopping. 



and we have introduced the Fourier transform of 

0(AA') ^ 



/^(A)|^(A') 



(10) 

(11) 



All these matrix elements can be expressed in terms of the 
coefficients aj"^^ as discussed in the Appendix VK\ Once 
these matrix elements are calculated we readily obtain 
the band structure for the spin polarons. 

To conclude this section we briefly discuss the relation- 
ship with previous work. Several authors have studied 
hole motion in an antifcrromagnet by calculations within 
a string basis 4, 45, 46, 47]. The difference is that the 
diagonalization of Hq leads to a considerable reduction 
of basis states in that high lying eigenstate of Hq are 
eliminated from the very beginning. The matrices to be 
diagonalized in the present work are 4x4 or lOx 10- 
which is very small compared to the matrix dimensions 
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in Refs. [1 EE [H, [13, Sg. Moreover the LCAO-like 
scheme makes it easier to extract a physical picture. An- 
other frequently appUed approach is the self-consistent 
Born approximation0, i, i, 111 [ul, [H H 113, [Hi • The 
single-hole wave function associated with this approxima- 
tion actually also can be interpreted as a superposition of 
string states once the Fourier-transformed version of the 
wave function in Ref. [52\ is converted into real space. 
This explains why the results e.g. for the dispersion of a 
single hole are practically identical. 
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IV. PHOTOEMISSION SPECTRA AND FERMI 
SURFACE 

ARPES gives the information on the one-electron re- 
moval part of the spectral function defined, at T — 0, 
as 

^ -Ck,J*Af). 

(12) 

I^^af) represents in the half-filled ground state in 
which the photoemission process takes place. Since we 
assume that the rigid band scenario is applicable to 
cuprates in the low doping range we expect that the con- 
clusions drawn from that analysis are also to some extent 
applicable to doped systems. 

As a first step we approximate the resolvent operator 
{{uj-H + iO+y^ by 



1 



l'i>;(k))(<i>i(k)| 



H + iO+ 



^0^ 



(13) 



i.e. we restrict the single hole states to the coherent su- 
perposition of SP states (O. 

Next, we have to choose an approximate ground state 
IvP^f) of the Heisenberg antiferromagnet. The simplest 
choice would be the Neel state \^n) but in this way we 
would miss an important mechanism for a k-dependent 
quasiparticle weight, namely the coupling of the photo- 
hole to quantum spin fluctuations. By generating a hole 
in the Neel state we can obtain only a bare hole, i.e. 
the string of length 0. In the presence of quantum spin 
fluctuations the photoemission process can also generate 
strings of length 1 or 2 - see Figure ID Since in such a 
process the hole is created not at the central site i of the 
SP state - which determines the phase factor e~^^'^' in 
the Bloch state (O - the photoemission matrix element 
becomes k-dependent. Generally speaking the fact that 
the SP quasiparticles extend over more than one unit cell 
in real space results in a 'structure factor' which varies 
within the first Brillouin zone in k-space. As will be seen 
below, this very k-dependence is the source of the rem- 
nant Fermi surface. In order to capture this effect we use 
simple first order perturbation theory for the quantum 
spin fluctuations and set 
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FIG. 4: (Color online) (a) (b): Creation of a bare hole at 
site i from the Neel state, (c) — > (d): Creation of a string of 
length 1 starting at i from the Neel state -I- quantum fluctu- 
ation, (e) — > (f): Creation of a string of length 2 starting at 
i from the Neel state + quantum fluctuation. 



Since second order perturbation theory gives quite a good 
estimate for the ground state energy of the Heisenberg 
antiferromagnet we expect that the probability for cou- 
pling to a quantum fluctuation with the electron annihi- 
lation operator is described quite well by ^T^ ■ 



(id) 



Since a quantum fluctuation in the initial states sim- 
ply gives rise to an extra factor of — 1/3 we immediately 
obtain the following expressions for the photoemission 
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matrix element rrS^'^ 



(*l'Hk)|ck„ 



i\^AF) 



TO(")(k) 



- COB{ky)) 



2ai , ^ 
ao (cos(fc^) + 

((*^°'^(^^) + cos(fcy)) 2-1) 

— ^ — sm(fea;) 

H g sm(fc3;) (cos(fc^) + cos{ky)) 

(cos(fc^) - cos(fcj,)) 

_^(cos2(fc,)-cos2(fc;,)) (15) 

Using these matrix elements we can now compute the 
spectral density from the normalized SP eigenfunctions. 

If we want to compare to experiment, however, there 
is yet another important effect we need to take into ac- 
count, namely the coupling of the photohole to charge 
fluctuations. We may expect that the ground state of the 
system has not only quantum spin fluctuations but also 
charge fluctuations i.e. an admixture of pairs of holes and 
double occupancies with a density oc {t/UY- By annihi- 
lating an electron on a doubly occupied site it is possible 
to create a string state with an initial site i by annihi- 
lating an electron at a site different from i. Again, this 
will give rise to a k-dependence of the spectral weight. 
Such processes actually are not described by the iJM, 
but since it is rather easy to discuss them we do so. We 
treat the charge fluctuations in perturbation theory i.e. 
we replace 



1^ 



AFI 



1^ 



AFI 



(16) 



where d, 



The hopping terms oc t',f" do not 
produce charge fluctuations in the Neel state. Denoting 
1] = t/U — J/{4:t) we get the following correction to the 
matrix element m^"^ 

m('')(k) = 277q;o (cos(fcj;) + cos(fcj^)) - 47701 (17) 

whereas the corrections to m^P^^ and m'^'^^ are zero. It 
might appear that the corrections due to charge fluctua- 
tions are quite small, being anyway oc t/U . On the other 
hand, by coupling to charge fluctuations a bare hole at 
site i can be created by actually annihilating an electron 
on any of its z neighbors. By changing from k — (0, 0) 
to k = (tt, tt) the corresponding contribution to the pho- 
toemission matrix element thus changes from zt/U to 
—zt/U and since this has to be added before squaring 
the matrix element the impact of the charge fluctuations 
is in fact quite strong. The way in which we are treating 
charge fluctuations would be adequate for a simple one- 
band Hubbard model, which is not the proper model for 
cuprate superconductors. A very similar calculation for 
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FIG. 5: Band structure obtained by solving the eigenvalue 
problem ([7]). Also shown is the disp ersion of bands observed 
by Ronning at al.m Ca2Cu02Cl2[54| (circles and triangles). A 
sample Fermi level of a doped system has been chosen as the 
zero of energy. Vertical arrows label possible optical transi- 
tions of a hole at the bottom of the band. 



the more correct two-band model has bee done by Eroles 
et o/.Hl]. 

In a first calculation we want to study the low energy 
band structure. In the LCAO-like ansatz ([6]) we first 
restrict ourselves to the lowest state (i.e. to = 1) for 
each symmetry o so that we have to solve 4x4 matrices 
(we have o € {s,Px,Py,dx'^_y-2^). The resulting band 
structure is shown in figure [S] Using the photoemis- 
sion matrix elements (jl5p the photoemission spectrum 
can be calculated, see figure [51 In addition to the fa- 
miliar quasiparticle band discovered in Sr2Cu02Cl2 by 
Wells at a/. [J, which has been discussed extensively in 
the literatureta, i, 0, i, i, E, [HI, [H, S HI there is a 
second band with slightly lower intensity which has pre- 
dominant p-like SP character and runs essentially paral- 
lel to the original band. Indeed such a second band with 
weaker intensity which follows the main band is seen in 
exact diagonalization of small clusters - see e.g. Figure 
1 of Ref 17]. These higher lying SP bands have small 
spectral weight and therefore could be hard to observe 
in ARPES - even more so because the ARPES spectra 
in the undoped compound are likely to show strong lat- 
tice polaronic effects [SSf. Nevertheless, it may be that 
a higher lying band - possibly the first p-like band - 
has been observed by Ronning et al. in the insulating 
cuprate Ca2Cu02Cl2[54'|. There a second weak band has 
been observed at roughly 0.5eV below the quasiparti- 
cle band. Ronning et al. interpreted this as part of a 
wide band which reaches F at a binding energy below 
2eV. As can be seen in Figure [51 however, this band 
(triangles) has a dispersion that is quite comparable to 
that of the quasiparticle band (circles). Moreover, as 
will be discussed below, we believe that the band portion 
observed in Ca2Cu02Cl2 around F at binding energies 
around 2eV is an example of a 'leV peak' as observed in 
Sr2Cu02Cl2[Il| and thus unrelated to the band marked 
by triangles. While the higher lying bands may be hard 
to observe in ARPES optical interband transitions be- 
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FIG. 6: Photoemission spectrum for the half-filled model cor- 
responding to the band structme Q. In the left panel the 
spectral weight is computed using only the quantum fluctua- 
tion correction (jlSp . in the right panel the charge-fluctuation 
correction ([TT)) was used as well. 
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FIG. 7: Top: Fermi surface obtained for 10% hole concentra- 
tion by rigid filling of the lowest SP-band. Bottom: Spectral 
weight of the corresponding SP-band along the Fermi surface 
as a function of the angle a. 



tween the main band and these higher lying bands pro- 
duce finite-frequency optical conductivity which may cor- 
respond to the mid IR bands - again in the actual com- 
pounds this may be complicated by polaronic effects. 
Another feature which can be seen in Figure [5] and which 
is quite consistent with experiment is the sharp drop of 
the spectral weight of the quasiparticle band which oc- 
curs whenever one passes (roughly) through the Fermi 
surface of the noninteracting half-filled, i.e. the rem- 
nant Fermi surface. It is caused by the k-dependence 
of the photoemission matrix elements and ([T7]) and 
therefore reflects the interplay of spin and charge fluc- 
tuations in the spin background and the internal struc- 
ture of the quasiparticle. Assuming that the structure 
of the quasiparticles remains roughly the same in a spin 
background without long range order but short ranged 
antifcrromagnetic correlations - as is suggested by ex- 
act diagonalization[32 - the k-dependence of the pho- 
toemission matrix element should be similar for flnite 
doping. This would provide an immediate explanation 
for the 'Fermi arcs' seen in the underdoped compounds. 
To make this more quantitative we have computed the 
Fermi contour by filling up the single hole dispersion ac- 
cording to the Pauli principle - sec Figure [5]- and showed 
it in the upper panel of Fig. [T] The lower panel shows the 
spectral weight of the quasiparticle band as a function of 
the Fermi surface angle. The hole pocket is actually more 
elongated along the (1,1) direction than along the antifcr- 
romagnetic zone boundary. It should be noted, however. 



that this is calculated at half-filling where (tt, 0) is far 
from the valence band top due to the t' and t" terms - 
see Figure [5] On the other hand it is known that upon 
doping the band portions near (7r,0) move upward(l3[. 
It has been suggested that decreasing antifcrromagnetic 
spin correlations lead to an effective downward renor- 
malization of t' and t" with doping [13] so that for finite 
doping the pocket is probably elongated more along the 
antifcrromagnetic zone boundary as observed by Chang 
et aL[20|. This effect cannot be reproduced by our sim- 
ple theory, however. Another phenomenon which would 
find a very simple explanation in an approximate rigid 
band behaviour upon doping is the pseudogapp^. The 
upper part of Figure [5] shows again the hole pocket and a 
contour in k-space which extends the 'front part' of the 
pocket into a free-electron-like Fermi surface. The pseu- 
dogap is usually defined by measuring the leading edge 
shift or the dispersion of the high-energy feature along 
such a contour. The lower part of Figure [5] then shows 
the energy of the quasiparticle band along this 'Fermi 
surface' plotted versus the Fermi surface angle <&. There 
is the characteristic flat part near = which originates 
because the contour initially follows the hole pocket and 
then the d-wave-like downward dispersion as the free- 
electron-like Fermi surface departs from the hole pocket. 
It should be noted that this would explain only the high- 
energy pseudogap. On the other hand the low-energy 
pseudogap, being defined in terms of a leading-edge shift, 
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FIG. 8: Top: Fermi surface obtained for 10% hole concen- 
tration by rigid filling of the lowest SP-band and a contour 
obtained by extending the 'inner part' of the pocket to a 
free-electron-like Fermi surface. Bottom: Energy relative to 
Fermi energy of the SP band forming the pocket along the 
free-electron-like Fermi surface as a function of Fermi surface 
angle. 



has no immediate connection with a dispersion relation 
and will almost certainly depend also on the tempera- 
ture and momentum dependent linewidth of the quasi- 
particle band, 18]. In fact, assuming a lifetime broaden- 
ing r(k) (X |i?(k) — /i| would immediately explain also 
the low-energy pseudogap. As noted above the pocket in 
Figure [5] is too much elongated in (1, l)-direction. How- 
vere, the effective downward renormalization of t' and t" 
in the doped compounds would lead to a pocket that is 
more elongated along the magnetic zone boundary. In 
any way, however, if one would go along the inner part of 
the pocket near (1,1) and extend this to a free-electron- 
like Fermi surface as in the top part of Figure [8] one will 
always see a dispersion as shown in the bottom part of 
Figure E] 

Finally we want to discuss the spectral function on 
a larger energy scale. To keep the discussion simple we 
keep only s-like SP states in ([6]) but include excited s-like 
states with m — 1 ... 10 (extending this to m = 1 ... 20 
produces no visible change in the spectra - this is one of 
the beneficial effects of the 'prediagonalization' of Hq). 
Moreover we retain only the matrix elements due to 
string truncation and the t' and t"-terms, i.e. processes 



of the type shown in Figure [3l However, we do include 
the full photoemission matrix elements (fTS]) and (|17p . 
The resulting spectral function is shown in Figure [9l In 
addition to the quasiparticle band - which is essentially 
identical to the one of the more exact calculation shown 
in Figure [6] above - there now appear additional bands 
at higher energy. In reality these 'bands' probably are 
not well-defined states because they are already high in 
energy. Rather these states are probably strongly broad- 
ened due to interaction with magnons and phonons and 
may have only the character of 'resonances'. Whereas 
the quasiparticle band is composed mainly of the low- 
est s-like SP states for motion of the hole trapped in 
the linearly ascending potential the higher lying bands 
correspond to excited levels of the trapped hole. The rel- 
atively high intensity of these states may be understood 
by noting that the coefficients ai™'' for the excited states 
m > 1 will have extra nodes as functions of v and if the 
signs of the alj^'^ better match the prefactors in the 
matrix element may even be larger for these higher lying 
states. 

The spectral function in Figure [9] is qualitatively similar 
to the result of a recent calculation by Bonca et aL\4T\ 
which was performed in a string basis with several mil- 
lion basis states. As already noted the main difference 
between the present calculation and the one by Bonca 
et al. is in the fact that the 'prediagonalization' of Hq 
by solution of ^ leads to a quite massive reduction of 
irrelevant degrees of freedom in the present scheme - as 
noted above, the matrix diagonalized here is 10 x 10. 
Despite this simplification not only the dispersion of the 
topmost peak but also the fact that the high energy part 
'widens' as one moves away from (0, 0) is reproduced. On 
the other hand, it should be also noted that the spectral 
weight shown in Figure [9] cannot directly be compared 
with Ref. [131 because we are taking into account the 
corrections of the ARPES matrix element due to quan- 
tum spin fluctuations in the half-filled gound state which 
are not included there. 



V. ELECTRONIC STRUCTURE AT HIGHER 
BINDING ENERGIES 

Recently a number of studies have revealed additional 
structure in the ARPES spectra at higher binding en- 
ergies. Ronning at al. found (5 4j that in the insulating 
cuprate Ca2Cu02Cl2 the quasiparticle band seem to 'fade 
away' as the F-point is approached - a phenomenon which 
is common to all cuprates. The spectral weight that 
is missing from the quasiparticle band then appears at 
~ 1.5eV below the maximum of the quasiparticle band 
in the form a high intensity band which has a dispersion 
that is remarkably consistent with LDA band structure 
calculations. More precisely this is the antibonding band 
of Cu3dx2-y2 and 02pa orbitals - whereby it has to be 
kept in mind that at F these two orbitals do not hybridize 
due to partity, so that in an LCAO-like description these 
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FIG. 9; Spectral density for hole creation at half-filling along 
high-symmetry lines in the Brillouin zone. Calculated with 
(right) and without (left) the contribution from charge fluc- 
tuations. 



bands have pure oxygen character at F and - by continu- 
ity in k - very small CuZd admixture its neighborhood. 

Similar behaviour was observed in doped cuprates 
as welljH^, m, [g^, [6l|. In addition, the so-called 
waterfall-phenomenon is observed. Moving e.g. along 
the (1,1) direction towards F the quasiparticle band first 
disperses away from the Fermi edge but then - at a mo- 
mentum of approximately (f , f ) - seems to bend down 
sharply and drop almost vertically down to an energy of 
~ leV below the Fermi edge. The apparent vertical part 
of the dispersion - the 'waterfalls' - can be seen as a hump 
in momentum distribution curves at binding energies in 
the range of Q.'beV l.OeV Upon reaching « leV be- 
low Ep the vertical parts then merge with two LDA-like 
bands of high intensity. 

Similar behaviour - namely small band portions with 
a free electron-like dispersion and high spectral weight 
near high symmetry points of the Brillouin zone - has 
been observed previously in the insulating compound 
Sr2Cu02Cl2 by Pothuizen at aL.SG]. The interpretation 
given by these authors was that these are 02p derived 
states which do not hybridize with Cu3d orbitals due 
to symmetry - which is why they appear only at high- 
symmetry points of the Brillouin zone - and hence are un- 
affected by the strong correlations in the partially filled 
Cu3d orbitals. This explains their LDA-like dispersion 
and high spectral weight because they are essentially free 
electron states. 

We believe that an important clue to the interpreta- 
tion of the waterfalls is the finding of Inosov et al. who 



showed that matrix element effects play a crucial role 
in their observation ,62J and that the missing part of the 
quasiparticle band near F can in fact be observed with 
photon energies around lOOey where the cross section for 
Cuid orbitals becomes appreciable , 63]. Moreover, Pan at 
al. found that the waterfalls show the same dependence 
on photon polarization as the quasiparticle band itself, 
indicating that they also are derived from Zhang-Rice 
singlets [60] . And finally in our opinion the crucial clue 
is the fact, that the quasiparticle band itself cannot be 
observed near the F-point either. While numerical stud- 
ies of the Hubbard and t-J modeljll, [H, [6^ - as well 
as the present theory - do indeed predict that the spec- 
tral weight of this band is lower by a factor of 2-3 at F 
as compared to (^, ^), in experiment there is practically 
no more intensity visible. Rather, the spectra show a 
complete suppression of spectral weight around F which 
extends down to the intense LDA-like bands. Most sig- 
nificantly, however, the waterfalls appear at very nearly 
the same momentum where the quasiparticle band itself 
becomes visible. 

We conclude that the reason for the vanishing of the 
quasiparticle band near F is 'extrinsic' to the t-J or single- 
band Hubbard model namely the special combination of 
phases for the 02pa orbitals in the bonding combination 
which hybridizes with a given CuSdx^_y2 orbital in the 
Zhang- Rice singlet (ZRS). In the framework of a simple 
three-step model of photoemission for photoemission the 
effect that we seek comes from the matrix element for 
the dipole transition from 02p states into the final state 
- which we take to be a plane wave with momentum k 
for simplicity. For definiteness we introduce the dipole 
matrix element 



1 



dr e 



r ■0a(r) 



(18) 



where e denotes the polarization vector of the light, V 
the volume of the crystal, V'a(r) is the wave function 
of an 02p orbital at the origin and a S {x,y}. The 
matrix element for a dipole transition from the bonding 
combination of 02pcr oxygen orbitals around a given Cu 
site j 



P. 



1 



into the plane wave state is 



P 



-P 



mzRs 



(19) 



(20) 

Note that the k dependence of the expression in square 
brackets comes solely from the interplay between the 
phase factors e^^ '' on the four oxygen neighbors of atom 
j and the relative phases of the orbitals in Namely 
any two oxygen orbitals whose position in real space dif- 
fer by one lattice spacing have a relative phase of (—1) in 
the ZRS - which would correspond to momentum (7r,7r). 
This k dependence is therefore completely independent 
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of details in the computation of the matrix elements Vq 
and therefore in particular independent of the photon 
polarization. Moreover it would stay the same if a more 
realistic final state wave function were used as long as 
this is a Bloch state with momentum k. It now can be 
seen that mzRs ^ as k ^ (2n7r, 2TO7r). Unlike the 
argument based on the dipole selection rule which was 
discussed by Ronning et aZ.^] the expression thus 
explains why the quasiparticle band has vanishing inten- 
sity not only at T but at all equivalent points in higher 
Brillouin zones as well. 

The above considerations apply only to hole creation 
on oxygen. For photon energies around 20eV is is well 
known, however, that the cross section for hole creation 
on oxygen is considerably larger than for hole creation in 
transition metal 3d orbitals[63|. Moreover, starting from 
a state one would reach a state. The latter is high in 
energy and hence has small weight in the ZRS and more- 
over the resulting spectral weight would be spread over 
several eV due to the multiplet splitting of so that the 
corresponding matrix element certainly is small. A d^ 
final state could be produced by photoemission from d^^ 
state but this state also is high in energy and therefore 
has small weight in the ground state. In any way the con- 
tribution of hole creation on Cu is strongly suppressed. 
For larger photon energies - around lOOey - Inosov at 
al. could indeed resolve the quasiparticle bandfl^l but 
only at (27r, 0). Obviously the larger photoemission cross 
section and the non-validity of the dipole selection rule 
at this momentum make it possible to see intensity from 
photo holes in Cu3dx2-y2 orbitals. Clear evidence that 
the intensity seen there is due to photohole creation on 
Cu is also provided by the strong variation of intensity 
around the Cu3p Cu3d threshold which is being used 
routinely to identify transition metal 3d states in other 
transition metal oxides[67]. 

One might then ask, where the spectral weight cor- 
responding to hole creation on oxygen goes at F. The 
answer is that this spectral weight is concentrated in the 
LDA-like bands which are observed around the F-point 
at approximately 1.5eV below the Fermi energy. Due to 
parity the 'antibonding band' of Cuid^^-y-i and 02pa 
orbitals actually has pure oxygen character at F. If a 
Bloch state of 02p orbitals is created with momentum 
(0, 0) this cannot couple to the ZRS singlet - see (PD|) - 
but has overlap 1 with these bandlike states. This is the 
same reasoning as given by Pothuizen et al. j56i] for the 
'leV peaks' in Sr2Cu02Cl2. The spectral weight transfer 
from the quasiparticle band to the 'leV peaks' at F as 
observed by Ronning et al. thus is strong evidence for 
the ZRS-character of the quasiparticle band. 

These leaves the question as to what is the interpreta- 
tion of the waterfalls. The above considerations suggest 
that these are simply the higher lying SP bands seen in 
Figure[9l These bands are likely to be strongly broadened 
due to interaction with spin excitations and possible also 
phonons although one of them may have been resolved 
by Ronning et aLfs^j (see Figure [S]). Near F the spec- 
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FIG. 10: Interpretation of the waterfall phenomenon. Around 
r photoholes on oxygen do not couple to the ZRS whence all 
t-J derived states have no spectral weight. Photoholes cre- 
ated with this momentum on oxygen instead propagate as 
(nearly) free-electron states which do not hybridize with the 
correlated Cu d-orbitals. Moving away from F the t-J bands 
become visible but lose spectral weight as well (dashed po- 
tions). The nearly free-electron states cease to exist because 
they now have appreciable hybridization with the correlated 
Cu d-orbitals. The higher lying SP bands thus are visible only 
in a small window in k-space i.e. the waterfalls. 



tral weight of these bands disappears because they are 
also 't-J-derived' and the matrix element for creation of 
a ZRS vanishes see Figure [TUl As one moves away from 
F the matrix element for creation of a ZRS increases but 
- as can be seen in Figure [9]- the spectral weight of these 
excited SP bands now quickly decreases. It follows that 
the spectral weight of these bands must go through a 
maximum as one moves away from F and this is our in- 
terpretation of the 'waterfalls': a number of essentially 
incoherent states which have a 'window of visibility' in a 
narrow range of momenta around (f , f ) and this window 
of visibility is seen as the hump in the momentum distri- 
bution curves. Since the suppression of these states near 
F is governed by the same matrix element of the ZRS as 
the quasiparticle band itself, it is moreover clear that this 
window of visibility 'opens' precisely in the same range 
of k where the quasiparticle band itself becomes visible - 
hence the apparent downward bending and the waterfall- 
like appearance of the spectra. 

A more quantitative description of this phenomenon 
obviously would have to start out from a three-band 
model so as to describe both, the coupling of a photohole 
to a ZRS and the existence of the nearly free electron 
states at high-symmetry points. This is out of the scope 
of the present paper and we therefore make no attempt 
for a quantitative discussion. There have been a number 
of attempts to explain the waterfall phenomenon within 
the t-J modeljill, [6l|. However, as was already noted. 
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the experimental spectra show a complete suppression of 
spectral weight around T which extends downward all 
the way to the intense LDA-like bands. This behaviour 
is not reproduced by t-J model calculations which show 
spectral weight - corresponding to the higher lying bands 
in Figure [9]- at too low binding energy. Moreover, the 
LDA-like dispersion of the high-intensity parts near F is 
not really reproduced 
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VI. OPTICAL CONDUCTIVITY 

Finally we turn to a discussion of the optical conductiv- 
ity thereby using the results of the simplified calculation 
which took into account only the lowest (m =1) state 
for each symmetry o. As is the case for atomic wave 
functions the s and c?-like SP states on one hand and the 
p-like SP states on the other have opposite parity and 
hence can have nonvanishing matrix elements of the cur- 
rent operator in between them. If we assume again that 
the quasiparticle band is filled with holes upon doping - 
which would occupy momenta around (-j, -7) - we thus 
expect optical interband transitions (as indicated in Fig- 
ure O which should be observable in the finite frequency 
optical response. This is defined as 

C7siq^0,u;) - V- |(^„b5(q = 0)|Vo)P 

^ — ^ (n 



Siuj - [E^ - £;o))- 



where \ipn) {En ) denotes the n-th eigenstate (eigenen- 
ergy) of the system (in particular, n = denotes the 
ground state). Also, jg, with 6 = x,y denotes a compo- 
nent of the current operator. 



j(q) = *^i„„e''i-(^'"+^")/2 [R™-R„ ]cl 



(21) 



Assuming a filling of the quasiparticle band we approxi- 
mate this by 



asiu;) = 2^^1|($,(k)b,|<I>i(k))|V 

LU 

k 1 = 1 



(22) 



where denotes the hole occupation of the quasiparticle 
(1 = 1) band. To evaluate this we need matrix elements of 
the current operator between the localized SP states 
{^^lisl'i'i') (with A = (o, to)). The main contribution 
to this matrix element are shown in Figure [TT] Starting 
from a bare hole - the 'string of length 0' - the current 
operator jy generates two strings of length 1 but with 
opposite sign. These two string states thus have precisely 
the right sign to couple to a Py-like SP state, see the sign 
convention in Figure O From this and similar processes 
we obtain the matrix element 



(s),0\ 



/j=0 



^aX+i). (23) 



FIG. 11: Application of the current operator jy to a bare hole 
in the Neel state (a) creates two strings of length 1 (b) and 
(c) with opposite signs. 
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FIG. 12: Optical conductivity at 10 % hole concentration. 



All matrix elements generated in a similar way and in- 
volving SPs of different symmetries or hopping at longer 
distances have been listed in the Appendix |B] The re- 
sulting optical conductivity at a hole concentration of 10 
% is shown in Figure [121 While the higher lying bands 
with p-character have small spectral weight and hence are 
difficult to observe in ARPES, they dominate the optical 
conductivity. We note that the present interpretation of 
the optical conductivity is consistent with that of Ref. 
[4^ and [111. One may expect that the excited levels 
of the trapped hole - or 'internal degrees of freedom' of 
the SP quasiparticles - are also seen in the d yna mical 
density correlation function. Vojta and Becker |34l| have 
calculated the dynamical density correlation function in 
a string framework similar to the one used by Bonce at 
al. and obtained convincing agreement with exact diag- 
onalization results. Within our approach one can natu- 
rally explain the distribution and spreading of the spec- 
tral intensity over a wide energy range, which has been 
observed in experiments 69, 70]. Due to the Brillouin 
folding in the AF state the center of the hole pocket at 
(7r/2, 7r/2) is a high symmetry point, s, d, andp SP states 
do not mix with each other exactly at this point. Thus, 
the transitions with similar intensity occur in the weakly 
doped system between the lowest predominantly s-like 
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and the first and the third predominantly p-like bands. 
Wider spreading of the high intensity region at higher en- 
ergies can be attributed to the admixture of p-wave SPs 
to the lowest band and transition to the predominantly 
d-wave second excited band. 



VII. DISCUSSION AND CONCLUSIONS 

In summary we have presented a theory for spin 
polaron-like quasiparticles. The basic idea is that a hole 
in an antiferromagnetically ordered 'spin background' is 
self-trapped which actually requires only short-range an- 
tiferromagnetic order. This will lead to a hierarchy of lo- 
calized states which may also realize different irreducible 
representations of Civ Since the Heisenberg exchange 
and the t' and t" terms have matrix elements between 
such self-trapped states on neighboring sites an LCAO- 
like description emerges where the role of the atomic 
or Wannier functions is played by the levels of the self- 
trapped hole. This leads to a multi band structure for 
the doped holes with the lowest of these being the famil- 
iar quasiparticle band observed in the insulating com- 
pounds and discussed extensively in the literature. Here 
we take the point of view that the simplest description 
for the underdoped compounds is holes being filled into 
this quasiparticle band. The fact that these self-trapped 
states extend over several unit cells in real space necessar- 
ily implies that they have an ARPES form factor which 
varies within the first Brillouin zone - hence the strong 
variation of the photoemission intensity of the quasipar- 
ticle band as a function of k which explains the rem- 
nant Fermi surface and the Fermi arcs seen in ARPES. 
Moreover, the pseudogap becomes a triviality within this 
picture. One of the higher bands of the effective LCAO- 
Hamiltonian may have been observed in ARPES in the 
insulator Ca2Cu02Cl2 and optical transitions between 
the resulting bands may explain the mid-infrared band 
in optical spectroscopy. There are probably complica- 
tions due to lattice polaron effects but here we neglect 
these although the simplicity of the present calculation - 
which never need matrices to more than 10 x 10 - would 
certainly allow to treat such effects as well. 

The main drawback of the present theory is the use of 
a spin background with antiferromagnetic order - which 
clearly is not realized in doped materials of interest. On 
the other hand all processes by which the hole propagates 
involve only spins in its immediate neighborhood. One 
may therefore expect that very similar processes would 
occur in a spin background with only short range anti- 
ferromagnetic correlations so that much of the present 
theory should apply in this case as well. The most im- 
portant effect we are missing with the present calculation 
is the closing of the pseudogap with both increasing tem- 
perature and increasing doping. Since this closing implies 
that the dispersion actually approaches that for the 'pure' 
t-J model - i.e. without t' and t" terms - this effect could 
be described by an effective downward renormalization of 



the t' and t" terms. The mechanism may be the decrease 
of spin correlations: in the Neel state the t' and t" terms 
can transport a hole 'completely coherent' i.e. without 
creating a spin excitation. As the spin correlation length 
decreases and reaches the 'range' of these terms in real 
space - two lattice spacings - the t' and t" will increas- 
ingly generate spin excitations when transporting a hole 
so that they change their net effect from coherent hole 
transport to 'excitation generating' hole transport. 
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APPENDIX A: SPIN POLARON MODEL 

We first derive the equations for the coefficients 
in (HI), thereby assuming We denote by 

the sum of all string states with v inverted spins, 
each multiplied by the proper phase factors according to 
figured] On a Bethe lattice the number of such states is 
= /(°) z{z - ly-^ where /(°) = 1 for o = s,d^2_y2 
and /(°) = 1/2 for o = Px,Py Our normalized basis 

— 1/2 

states are i^) — Til/ |<i>^) which obey 

V 

because each of the n^+i basis string states in the bra is 
generated exactly once from a state in the ket and the 
matrix element for hopping of a hole is positive. The 
states ([2]) can be written as 

with jiy = rh!"^ OLh''^\ Performing the variational proce- 
dure - thereby using ([T]) - we obtain the set of equations 

Eof3o + V^tPi = EPo (Al) 
Ei(3i + ^/lt(3o + VI^t(32 = E(3i (A2) 
E^/3^ + Vz-U{l3^-i+P^+i) = Ef3^ (A3) 

For o — Px,Py, dx2-y2 one has to set Po — and discard 
the first equation. After introducing a cutoff for v these 
equations can be solved numerically. 

Now we proceed to analyse with some sample processes 
which give rise to nonvanishing matrix elements (jlOp and 
PT|) . The hopping integral t is the highest energy scale 
among all model parameters t, t', t" , and J. On the other 
hand, the exchange energy grows fast with the number 
of fluctuations. That number is directly related to the 
length of paths Vi . This mechanism is responsible for the 
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tendency towards hole confinement and the construction 
of SPs also relies on it. Quasiparticle deconfinement oc- 
curs due to processes which are mediated by hopping to 
second and third NN and by the action of the XY term 
in the exchange interaction. Only string states \'Pi) with 
a low number fluctuations can be involved in those pro- 
cesses at the low energy scale. Thus, in order to find ma- 
trix elements (jlOp we need to determine coupling between 
short-string states [Pi). That couphng is induced by the 
perturbation part Hi which by definition contains terms 
related to hopping to second and third NN sites and 
which also contains the XY term. We restrict our con- 
siderations to processes involving string states [Pi) which 
are related to paths not longer than 2 lattice spacings for 
s-wave SPs and to paths not longer than 3 lattice spac- 
ings for SPs with lower symmetry. In the latter case we 
take into account longer strings because the zero-length 
string state \Vi) does not contribute to the wavefunction 
([2|) of SPs with lower symmetry and thus the weight is 
shifted to states representing longer paths. The abso- 
lute values of prefactors , a'^ corresponding to lowest 
eigenstates of (|Aip - (|A3P rapidly decrease with the string 
length /i which additionally justifies the restriction of our 
considerations to strings with short length. At the level 
of the approximation which we apply, there exist more 
than 20 categories of processes which contribute to (fTO|l 
and (|lip . The differences between different categories 
concern the underlying mechanisms or the geometry of 
involved strings. Since the mechanisms which give rise 
to coupling between s-wave SPs were discussed in detail 
in the past ji^, ISl here we will mainly concentrate on 
the issue how the lowering of SP symmetry influences the 
coupling between SPs and we will discuss representative 
examples of processes which give rise to the hybridization 
between SPs. In order to keep the Hilbert space as small 
as possible we consider at first the low energy sector and 
analyze only lowest SP states with given symmetry. 

For example. Fig. [13] (a), (b) depicts a process during 
which a SP polaron is shifted to a second NN site, (a) 
represents a string state [Pi) which is a component of the 
wavefunction ([2]) of a SP created at site i. This state has 
been obtained by creating a bare hole in the Neel state. 
It appears with the prefactor ag in the superposition ([2]) 
defining the s-wave SP at site i. With the same prefac- 
tor appears the string state depicted in Fig. [13] (b). It 
is a component of the s-wave SP created at site j. Since 
string states (a) and (b) are coupled by the second NN 
hopping term a non-vanishing matrix element (jlOp be- 
tween s-wave SPs is generated. The related contribution 
to that matrix element is given by 

^^((..o)(..o))^^^^^^^^i)^,,^2, (A4) 

Within a convention for presenting matrix elements this 
contribution can be written as 

ss,x-|-y,r,W : d. (A5) 

The presence of the group symbol C4 in ()A5P means that 
additional contributions to pop can obtained by apply- 
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FIG. 13: (Color online) (a) and (b) A process which gives rise 
to the shift of the s-wave SP to a second NN site, (c) and 
(d) A process which generates corrections to eigenenergies E\ 
and E[ of SPs. 

ing to (|A4[) elements of C4 different than identity. Since 
ckg = states depicted in Fig. [13] (a) and (b) do not con- 
tribute to SP wavefunctions ([2]) with lower symmetry. 
Thus the above discussed lowest-order process does not 
generate hybridization in which such SPs are involved. 
By process order we understand, in terms of the conven- 
tional perturbation theory, the total number of Hamilto- 
nian actions necessary to transform a state representing 
a bare hole created in the AF background into another 
such state. It is clear that this parameter is related to 
the sum of lengths for paths involved in a given process. 

The origin of the formula for the integral ril' in (jA4p 
is rather obvious. It is the product of prefactors which 
appear by string states depicted in Fig. [13] (a) and (b) 
in the definition for the wavefunctions of SPs at sites i 
and j. t' is the integral which appears as a prefactor in 
the hopping term coupling those string states. The same 
scheme can be applied to deduce the form of hopping 
integrals which appear in further terms in 7^°° "* . 

Fig. [13] (c) and (d) depicts a process which gives rise 
to corrections to SP eigenenergies Ei and E[ of SPs 
and amendments to the diagonal term in ((SJ). Those 
corrections originate with the coupling between differ- 
ent string states generated by the hopping to second NN 
sites. Those string states contribute to the wavefunctions 
of SPs created at the site i. The paths Vi correspond- 
ing to string states Vi depicted in Fig. [13] (b), (c) have 
been shown in Fig. [TJ] (a), (b), respectively. From now 
on we will also use the latter form of diagrams to rep- 
resent string states. The same mechanism as discussed 
above works for pairs of longer strings which have identi- 
cal forms with the exception of hole positions at ends of 
them. The corrections to energies of s-wave d-wave and 
p-wave SPs are ril'^'', T^f^ and t^'^ respectively. Their 
explicit form can be found below. Due to point-group 
properties the above discussed process does not generate 
the coupling between SPs with different symmetries. 
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FIG. 14: (a) and (b) A different representation of string states 
depicted in Fig. [13] (c) and (d) respectively, (c) A string state 
contributing to a process which involves hopping to second 
NN sites. That process gives rise to SP shifts along plaquette 
diagonals, (d) A string state contributing to a process which 
involves hopping to NN sites. 

During another process which also involves strings of 
length 2 lattice spacings the action of hopping to second 
NN sites generates shifts of SPs between second NN sites. 
Surprisingly, SPs are shifted now in the direction opposite 
to the direction of the bare hole move. The nature of that 
paradox can be understood by analyzing Figures [14] (b) 
and (c). The action of the second NN hoping term on 
the string state, Fig. [M] (b), which contributes to the SP 
created at site i, shifts the hole at the end of that string 
from site j to site i and generates a string state. Fig. [14] 
(c), which contributes to the SP created at site j. That 

process gives rise to a new contribution to , 

ss, X + y, T^l^^ ; dd,x + y, r^f^ ; 

Pi;Px,x + y,r^p^' : C4V (A6) 

Pys, i + y, T^i5) ; pyd, X + y, r^^'^ : C4. • H.C{A7) 

The presence of the group symbol Ciy means as before 
that the contributions (jA6|) and (jA7|) should be supple- 
mented by matrix elements obtained by applying to the 
existing elements all group elements and in the case of the 
second line also elements obtained by applying the Her- 
mitian conjugation. The process depicted in Fig. [TJ] (b) 
and (c) gives rise to the hybridization of SPs with differ- 
ent symmetries because apart from the case of coupling 
between s-wave and d-wave SPs no selection rule forbids 
it and also because string states depicted in both Fig- 
ures [14] (b) and (c) can be components of wavefunctions 
for SPs of arbitrary symmetry. The hopping integrals 

T^Q, as usually are given by products of the bare hop- 
ping integral t' and of prefactors with which string states 
depicted in Fig. [TJ] (b) and (c) appear in the definition 
of SPs. The existence of negative signs in some of inte- 

(15) 

grals T^Q, originates with the fact that due to symmetry 



properties of some SPs such signs are also contained in 
prefactors of string states in sums ^ defining those SPs. 

After discussing some representative examples of SP 
coupling generated by hopping to second NN sites we 
skip the discussion of analogous processes mediated by 
hopping to third NN sites, because their mechanisms are 
similar, and analyze a high-order process which involves 
hopping to NN sites. That process has been neglected 
during the construction SPs, because a long-string state 
participates in it. The NN hopping operator, by shifting 
the hole at the end of the string depicted in Fig. [TJ] (d) 
from site n to site j, generates the string state depicted in 
Fig. [13] (b), which gives rise to a new contribution to the 
coupling between SPs created at sites i and j. For the 
sake of clarity, we only additionally mention briefly, that 
Fig. [14] (d) represents a string state in which flipped spins 
occupy sites j, I, i, and the hole occupy site n. It is clear 
that the same mechanism generates the hopping of SPs 
between all second NN sites. Depending on the symmetry 
of the initial and final SP states, related hopping integrals 
are given by T^f\ T^f\ Tj^f\ t^I^\ and 4f\ 

The analysis of Hamiltonian terms induced by the ex- 
change interaction can be performed along similar lines. 
We start with a necessary correction to the Ising part of 
the exchange energy. When we have been constructing 
the Schrodinger equation (|A1[) - (|A3|) for wavefunctions ([2]) 
of SPs, we have assumed that the Ising contribution de- 
pends on the length of strings, but does not depend on 
their geometry. That assumption is true for short strings, 
but for longer strings there exists some exceptions from 
that rule. For example, by inspecting the spin structure 
of the string state depicted Fig. [TJ] (d), we see that one 
broken bond has been saved for that state in compari- 
son with strings which have the same length and have 
the starting and ending points which are not NNs. That 
fact gives rise to the decrease by J/2 of the string-state 
energy and to corrections to eigenenergies of d-wave and 
p-wave SPs. Those corrections are given by and tpp^'' 
respectively. Since the amendments are related to strings 
of length at least 3 lattice spacings we neglect them in the 
case of the s-wave SP because the weight of related string 
states in the wavefunction of the s-wave SP is small. 

Also the best known in literature process which de- 
termines to great extent the overall shape of the single- 
particle energy-dispersion at the low energy scale gives 
rise to the hybridization between SPs with different sym- 
metries. That process has been shown in Fig. [15] (a) and 
(b). Flipping spins at sites I and m in the string state 
depicted in Fig. [T^] (a) by the transverse part of the ex- 
change energy gives rise to the string state depicted in 
Fig. [15] (b) and to the hybridization between SPs at sites 
i and m. In that kind of hybridization the s-wave state 
should be involved because the state shown in Fig. [TSlfa) 
represents a bare hole created in the Neel state. As we 
have shown before, such string states do not contribute to 
wavefunctions of SPs with lower symmetry. The hopping 
integrals which appear at terms coupling an s-wave SP 
with an s-wave, a c?-wave, and a p-wave SP located at a 
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FIG. 15: (a) and (b) A low-order process which involves the 
action of the transverse part of the exchange term and gives 
rise to the hybridization of SPs with different symmetry, (c) 
and (d) String states involved in a process mediated by the 
XY part of the exchange term and resulting in a correction 
to the eigenenergy of the d-wave SP. 



third NN site are given by t. 



(3) ,(3) 



, and i^s respectively. 



FIG. 16: (a) and (b) Spuriously different string states which 
are actually identical, what gives rise to the overlap between 
SPs created at sites n and I. 



gives rise to a new contribution in the Hamilton operator, 
Pyd, X + y, E[iu^^J : ■ H.C. (A8) 
dd, X + y, E[JJ-J ■ pyp, , X + y, E'^ujg^ : C4. . (A9) 

The analysis of other processes which give rise to new 
elements of ^ is rather straightforward and similar 
to the analysis of previously discussed processes. Thus 

we do not discuss the remaining contributions to ■* 
one by one. 

Finally, we list here matrix matrix elements which de- 
termine the form of the Hamilton matrix (HI), (flOl). 



A similar mechanism gives rise to a correction to the 
eigenenergy of a d-wave SP. The action of the XY term 
flips spins in the string state depicted in Fig. [15] (c) at 
sites n, j and creates the string state shown in Fig. [T5] 
(d). Due to the fact that both string states contribute 
to the wavefunctions ^ of SPs created at the same site 
i and due to standard selection rules, the contribution 
to '7^°° '' originating with the above discussed process 
must be diagonal. For this reason, p-wave states can 
not contribute to the related new matrix elements ' , 
because the process depicted in Fig. [15] (c) , (d) involves 
string states which contribute to different p2;-wave and 
Py-wave SPs. Since the length of the string depicted in 
Fig. [15] (d) is 3 lattice spacings we neglect the correction 
to the eigenenergy of the s-wave SP, because we expect 
that it is small. The correction to the eigenenergy of the 

d-wave SP is —i^^i^P ■ It also contains contributions related 
to the coupling, in the same way, between longer-string 
states obtained by letting holes at the ends of strings in 
Fig. [15] (c) and (d) to hop further along identical paths. 

During the construction of SPs we have been assuming 
that wavefunctions of SPs at different sites are orthogonal 
which turns out not to be an exactly true assumption. On 
the other hand the overlap between SPs is rather small 
because it originates with the overlap of nominally dif- 
ferent string states related with paths of relatively long 
length. Such string states have been depicted in Fig. [TC] 
(a) and (b). Both Figures [16] (a) and (b) actually repre- 
sent the same state which consists of a hole at site j and 
fluctuations at sites n, i, I created in the Neel state. The 
generated in this way overlap between different SPs also 
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and of the overlap matrix ([9]), ((TT|) . 



where 
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For example contributions No. 1 and 2 are related to 
shortening and elongating strings by the current opera- 
tor, Fig. [Til while contributions No. 3 and 4 to shifts 
between next NN sites, Fig. [rHTa) and (b). We have con- 
sidered processes involving strings of length up to two for 
matrix elements of the t-tcrm in the current operator and 
up to one for matrix elements of t' and t" terms. 



APPENDIX B: OPTICAL CONDUCTIVITY OF 
DOPED ANTIFERROMAGNETS 



The optical spectrum evaluated by us is deter- 
mined by following contributions to matrix elements 

'5^"^*r;*"|>I*r=o) = si>(R») where n labels differ- 
ent contributions, 

4'is(0) = -2zt^a,,a;+i/y2, 
4K0) = -2it^a;a^+i/\/2, 

(0) = -2»i^aX+i/(2V2)> 

si'yo) = s«JO) (Bl) 
■(2),(0) = 2zf^a^a;,_i/y2, 



V. 

s 



41(0) = 2ji^aX-i/^> 

41(0)--2*tE"X-i/(2\^)' 

4!!<^(0)=41(0) (B2) 

43)(x±y) =-zt'a2 (B3) 

si%~i.±y) = ~si%i + y) (B4) 
45.(0) ==4i^Vlal/^/2, 

41(0) = ~4*i'K)V(2V2) (B5) 

4S(0) = 41(0) (B6) 
4yi(2x) = -2zi"a2, 

.S(-2x) = -.S(2x) (B7) 

Sp],sW = -4:it"aia[/V2, 

41(0) = -41(0), 

41(0)=4zt"K)V(2V2), 

41(0) = -41(0). (B8) 
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